Mass ejection by strange star mergers and observational implications 
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We determine the Galactic production rate of strangelets as a canonical input to calculations of 
the measurable cosmic ray flux of strangelets by performing simulations of strange star mergers and 
combining the results with recent estimates of stellar binary populations. We find that the flux 
depends sensitively on the bag constant of the MIT bag model of QCD and disappears for high 
values of the bag constant and thus more compact strange stars. In the latter case strange stars 
could coexist with ordinary neutron stars as they are not converted by the capture of cosmic ray 
strangelets. An unambiguous detection of an ordinary neutron star would then not rule out the 
strange matter hypothesis. 
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The strange matter hypothesis (SMH) considers the 
possibility that the absolute ground state of matter might 
not be formed by iron but by strange quark matter 
(SQM), a mixture of up, down and strange quarks P, 0. 
If true, stable objects consisting of this matter with 
baryon numbers from about 10^ to 10^^ might exist 
0j 0, Hi- The latter end corresponds to strange stars 
(SS) with a mass and radius comparable to that of neu- 
tron stars (NS), where the upper mass limit is given by 
the inevitable collapse to a black hole (BH) p, 0I- I^i 
contrast to NSs these SSs are selfbound and do not have 
an overall inverse mass-radius relation. 

One of the astrophysical consequences of the SMH is 
the possibility that collision events of two SSs lead to the 
ejection of strangelets, small lumps of SQM These 
strangelets would contribute to the cosmic ray flux. The 
pollution of the Galaxy with strangelets is speculated to 
convert all ordinary NSs to SSs. It was argued that all 
compact stars have become SSs in this scenario because 
already a tiny amount of strangelets is sufficient to trigger 
the transformation 0, [l^l ■ If this sequence of arguments 
was true, the unambiguous observation of a NS would 
rule out the SMH according to references Q . 

Estimates of the strangelet flux use results of NS-NS 
merger simulations [ll|, which arc not necessarily reli- 
able in the case of SQM. No detailed simulations of SS 
coalescence have been conducted so far, and the ejected 
mass is unknown. Only Newtonian simulations of SS-BH 
binaries, modeling the BH by a pseudo-relativistic po- 
tential [iJl , were carried out by [l3|. It was found that 
from this kind of mergers no matter is ejected. In or- 
der to shed light on the merger process of SS binaries we 
performed relativistic three-dimensional hydrodynamical 
simulations of the coalescence. 



Several current and upcoming experiments have the 
potential to detect signatures of SQM. For instance SQM 
might be produced directly in the Large Hadron Collider 
at CERN 
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15|. But also cosmic ray experiments like 
the Alpha Magnetic Spectrometer AMS-02 planned to 
be installed on the International Space Station in 2010 
are designed to capture strangelets [lil . 17 1. In addi- 
tion, gravitational- wave detectors like LIGO and VIRGO 
might identify characteristic signals from SS mergers and 
SS oscillations or instabilities Also indirectly, the 
observation of compact stars can help to decide on the 
SMH^s_^ecially by pinning down the mass-radius rela- 
For a review on additional SQM searches 

The expectation that all NSs convert to SSs and that 
there is a measurable flux of strangelets as cosmic rays, 
relies on the assumption that SS mergers or another 
source indeed eject SQM in a sufficient amount. Here 
we report that a pollution through SS mergers does not 
need not be present for all models describing absolutely 
stable SQM. In fact, we find that the amount of ejected 
matter depends on the so-called bag constant, which in 
the MIT bag model adopted here represents the pres- 
sure of the non-perturbative QCD vacuum. Therefore 
the determination of the mass flux of strangelets in cos- 
mic rays could help to constrain this unknown param- 
eter, which in turn gives the binding energy of SQM. 
Complementary insights in the equation of state (EoS) 
of SQM might come from the detection of gravitational 
waves (GWs) from SS coalescences, if some specific fea- 
tures distinguished them from those of merging NSs. 

In order to describe the EoS of absolutely stable SQM 
we employ the MIT bag model [H Within this 

model, quarks are considered as a free or weakly inter- 
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acting Fermi gas and the non-perturbative QCD inter- 
action is simulated by a finite pressure of tfie vacuum, 
the bag constant B. The small current masses of up and 
down quarks allow us to treat them as massless parti- 
cles, whereas for the strange quark we adopt the value 
of TOs = 100 MeV [2l|. For our study we consider free 
quarks, which corresponds to a range of the bag con- 
stant of 57 MeV/fm3 < B < 84 MeV/fm^. The lower 
limit of B is given by the fact that baryons do not con- 
vert spontaneously to a two flavor quark phase, while 
the upper limit is determined by the requirement of ab- 
solutely stable SQM (energy per baryon at zero pressure, 
E/A, lower than the corresponding value of 930 McV for 
nuclear matter). These limits can be altered for other 
choices of or by considering the interactions among 
quarks. The values of 60 MeV/fm^ {E/A = 860 MeV) 
and 80 MeV/fm^ {E/A = 921 MeV) for the bag con- 
stant have been chosen to represent the extreme cases of 
the underlying microphysical model, which we refer to as 
MIT60 and MIT80. This choice of parameters yields a 
maximum mass of bare cold SSs of 1.88 Mq for MIT60 
and 1.64 AIq for MIT80 with corresponding radii of 10.4 
km and 9.0 km, respectively. For given mass the stellar 
radii for MIT80 are in general slightly smaller than those 
for MIT60. A possible nuclear crust of SSs is neglected 
because of its small mass ('^lO~^A/0), which makes it 
irrelevant for the dynamics of the system. 

We performed SS merger simulations with the code de- 
scribed in [i^] . A three-dimensional relativistic smoothed 
particle hydrodynamics scheme (SPH) is combined with 
an approximate treatment of general relativity (GR) em- 
ploying the conformal flatness condition, supplemented 
by a method accounting for the backreaction of GW emis- 
sion on the fluid [24I. Magnetic fields are not included. 



because they can be considered as unimportant for the 
dynamical behaviour as long as the initial field strength 
inside the compact star is below '-^lO^^ G psj . 

The models of our simulations are chosen such that 
they cover the whole potential mass range of compact 
star binaries. The gravitational masses of the stars vary 
between 0.9 Mq and roughly the maximum SS mass for 
each EoS. The simulations start after a relaxation phase 
from a quasi-cquilibrium orbit about two revolutions be- 
fore the actual merger. SQM has a shear viscosity com- 
parable to that of nuclear matter [241 (and with color- 
superconducting phases even lower (25[.), which is ex- 
pected to be too low to yield tidally locked systems [2^; 
therefore we consider irrotational configurations. Ther- 
mal effects are also taken into account as these were 
shown to be important for merger simulations with nu- 
clear EoSs, in particular when the ejection of mass is of 
interest |22| . 

In total we discuss results of 29 simulations for MIT60 
and 19 for MIT80 with a resolution of about 130,000 SPH 
particles. Using nonuniform particle masses we achieve 
a mass resolution of roughly 10~^AfQ in these runs. The 
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FIG. 1: Computed models for MIT60 and MIT80 in the Mi- 
A'f2-plane of the gravitational masses of the SS binaries. Filled 
circles denote prompt collapse to a BH, while open circles 
indicate the formation of a HMO. 



results were tested for convergence by additional simula- 
tions with higher resolution. 

There arc two possible outcomes of SS mergers. For 
relatively high masses of the binary components the 
merged object collapses promptly to a BH shortly after 
the stars come in contact. The forthcoming formation of 
a BH is indicated by a steep decrease of the lapse func- 
tion. Also the central density increases within a fraction 
of the sonic timescalc to values of twice the maximum 
density of a single nonrotating SS and thus exceeds the 
maximum density of stable, uniformly rotating "super- 
massive" SSs. If the masses are lower, the merger rem- 
nant can be transiently supported against collapse mainly 
by differential rotation. Such a "hypermassive object" 
(HMO) [i^l emits a characteristic GW signal, which is 
sensitive to the total mass of the binary and the EoS (see 
[2^ for NS mergers). This will be analyzed in a separate 
publication. Since the system mass is much larger than 
the mass limit of supermassivc SSs for the given EoS, the 
remnant collapses to a BH after the angular momentum 
has been redistributed [i^. 

Figure [1] gives an overview of the simulated binary 
mass configurations and their outcome. Filled circles in- 
dicate prompt collapse to a BH while open circles corre- 
spond to the formation of a HMO. Mi and M2 refer to 
the gravitational masses of the SSs in isolation. 

To estimate the amount of matter that becomes grav- 
itationally unbound, we use the criterion defined in [23 |. 
It considers the energy of a fluid particle in a comoving 
frame and applies if pressure forces are small in com- 
parison to gravitational forces, which is well fulfilled for 
particles leaving the merger site. In addition we cross- 
checked these results by a simple criterion that monitors 
how much matter expands away from the merger site. 
The ejecta estimates agree within a factor of less than 
two. Ejecta from HMOs originate from the tips of tidal 
tails that develop on timescales longer than the timescale 
of prompt collapse to a BH. In the case of such a prompt 
collapse no angular momentum can be redistributed from 
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FIG. 2: Ejected mass per merger event, color-coded and mea- 
sured in Mq, as function of the mass ratio q — M1/M2 and 
the total system mass Mtot ~ Mi + M2 of the binary con- 
figurations for the MIT60 and MIT80 EoSs. The white line 
separates binary mergers with and without ejecta. 



the center to the outer parts of the merged object because 
the matter in the inner part is swallowed quickly by the 
BH 23 . Thus particles potentially forming an accretion 



torus around the BH have no chance to end up in tidal 
tails and to gain enough energy to become unbound. 

Figure [2] shows the estimated amount of unbound mat- 
ter for given mass ratios q = M1/M2 and total binary 
masses Mtot = M1+M2 computed for MIT60 and MIT80 
by means of our ejecta criterion. Above a certain Mtot 
value we cannot determine any amount of ejecta (see 
white lines in Fig. [2]). If Altot is below this limit, we 
obtain a steep rise of the ejecta mass in a narrow region 
of the Mtot-9-plane in both EoS cases for q < 0.85. For 
MIT60 the region where more than O.OIMq of matter be- 
come unbound is located around a total mass of about 2.5 
Mq. For MIT80 significantly lower total masses are re- 
quired to obtain unbound matter and the ejected masses 
are lower as well. This dependence on the bag constant 
originates from the fact that MIT80 leads to more com- 
pact stars with correspondingly smaller radii, which im- 
pedes the tidal disruption. 

Since the ejected mass is very low in comparison to 
the system mass, we found a dependence on the chosen 
resolution and the initial setup of the SPH particles. The 
values of the ejected mass are uncertain within a factor 
of ~2. However, our conclusion that some configura- 
tions do not eject matter relies on the occurrence of a 
prompt collapse to a BH. This is a safe result of our sim- 
ulations within the employed approximations. Therefore 
the border between systems that eject matter and those 
that do not can be considered as well determined (see 
Fig. [2]). Only for equal-mass binaries the borderline in- 
cludes configurations that do not collapse promptly and 
still do not eject matter, because such systems do not 
form pronounced tidal arms (see 1221 for NS mergers). 

Population synthesis studies [29| provide probability 
distributions of compact star binaries dependent on their 



system parameters (e.g. q and Mtot)- Folding our results 
for the ejecta masses with these probability distributions 
allows us to estimate the ejected mass per merger event 
averaged over the whole population. These numbers can 
be used to derive more accurately the expected flux of 



strangelets in a detector like AMS-02 Assuming 
that the results of [2^ hold also for SSs and not only for 
NSs, wc compute for MIT60 a population-averaged ejecta 
mass of 8 • 10~^Mq. The uncertainties due to the limited 
resolution and the criterion for determining ejecta masses 
can change this result up to a factor of ~4. For MIT80 
we do not find any ejecta because only configurations not 
present in the adopted population eject matter. 

For a rough assessment of the uncertainties associated 
with the theoretical population synthesis studies we em- 
ployed a second data set based on observations of massive 
progenitor stars in double systems [lOl ■ Using theoretical 
results for the mass relation between NSs and progenitor 
stars [sH and ignoring complications due to binary evolu- 
tion effects, we derive a probability distribution function 
of compact binaries. Taking into account uncertainties in 
the determination of stellar masses, we obtain an average 
ejecta mass per event of (1.4 ... 2.8) ■ lO"" Mq for MIT60 
and again a vanishing ejecta mass for MIT80. 

The bag constant B is the only parameter varied be- 
tween the EoSs and determines the mass-radius relation 
of SSs as the crucial property for the merger dynamics 
" 2^ . 3^ . For intermediate values of B we expect smaller 



ejecta masses than MIT60 but higher than MIT80. The 
borderline between models with and without ejecta would 
then be shifted to an intermediate location as well. 

QCD perturbative corrections can be absorbed in an 
effective bag constant that can be chosen to yield mass- 
radius relations which agree well with the bag models we 
used [s^l • Color superconductivity has only a small effect 
on the EoS 
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l35|. However, quark interactions ch ang e 

the window for absolutely stable SQM 0, 0, [H, [37 1 . 

Our findings have important observational implica- 
tions. The mass-radius relation (in our study determined 
by the bag constant) strongly affects the amount of mat- 
ter ejected from SS mergers. Therefore a measured mass 
flux of strangelets would constrain this relation if quark 
star mergers were the main source of strangelets. A rel- 
atively high flux would be an indicator for less compact 
SSs, while no or only a low flux would only be consis- 
tent with more compact SSs. This would also put lim- 
its on the bag constant and so the binding energy of 
SQM. Assuming a Galactic merger rate of SS binaries 
of 10~^...10~* yr~^ [2^, our population-averaged ejecta 
mass of ~1O~*M0 for MIT60 yields a Galactic strangelet 
production rate of M = 10^^...10^^ yr~^. Since the 
flux of strangelets near Earth depends linearly on M, we 
derive a 10 to 100 times larger value than in (ll| . 

Even more relevant are the consequences if there are 
no other production mechanisms of strangelets. Our re- 
suhs for MIT80 imply that the SMH cannot be ruled out 
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but would only be compatible with compact SSs, if ex- 
periments like AMS-02 could not find any evidence for a 
nonzero strangelet flux. A strangelet flux below a critical 
limit would mean that no or not all NSs might have con- 
verted to SSs by capturing a strange nugget. SSs might 
then still form by nucleation of SQM drops, e.g. during 
stellar core collapse and explosion or by mass accretion 
of NSs in binaries when the central conditions reach some 
critical threshold for the phase transition to quark mat- 
ter 0, [1, [13, IS] ■ In this scenario there is a limiting mass 
above which SSs are formed while NSs exist below (see 
also [1^ ) . Thus in the case of a large bag constant NSs 
and SSs could be in coexistence. In the light of our sim- 
ulations the unambiguous observation of a NS would not 
rule out the SMH contrary to the suggestion in [1,0]. 

We stress that these conclusions from our simulations 
hold only if SS mergers are the only efficient sources 
of strangelet ejection. In fact, several other suggestions 
have been made, e.g. core-collapse supernova explosions 
[iol or the ejection by electric fields from the surface of a 
SS 41 1 if SQM nuggets were embedded in the crust [i^ ]. 
Despite the remaining uncertainties of our simulations 
like the approximate treatment of GR, which also docs 
not allow us to follow the formation of the BH, the lim- 
ited mass resolution, the simplified EoS, and the omission 
of magnetic fields and a nuclear crust, we expect that a 
more sophisticated approach will only yield quantitative 
shifts, changing the exact values of the ejecta masses and 
possibly insignificantly moving the border between con- 
figurations with and without cjccta. 

Our results might also apply to other forms of self- 
bound matter like pion-condensed nucleon matter [^, [4^ 
provided the stellar properties are similar. 

Future investigations, which should preferably be done 
in full GR, should consider EoSs including quark interac- 
tions and color superconductivity, or should use descrip- 
tions beyond the MIT bag model. Also SS-BH mergers 
should be reexamined in the GR framework. Finally, GW 
signals from mergers besides the stellar cooling behaviour 
[4 [4^ may be promising means to decide on the SMH 
or other forms of selfbound matter, once wave measure- 
ments will become possible. 
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